A magnetic analog of the isotope effect in cuprates 
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We present extensive magnetic measurements of the (Caa;Lai-2,)(Bai.75_a;Lao.25+2^)Cu30j/ system 
with its four different famihes (x) having a r™''^(x) variation of 28% and minimal structural changes. 
For each family, we measured the Neel temperature, the anisotropies of the magnetic interactions, 
and the spin glass temperature. Our results exhibit a universal relation Tc = cJus for all families, 
where c ~ 1, J is the in plane Heisenberg exchange, and Us is the carrier density. This relates 
cuprate superconductivity to magnetism in the same sense that phonon mediated superconductivity 
is related to atomic mass. 
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The critical temperature for superconductivity Tf. in 
metallic superconductors varies with isotope substitution 
This observation, known as the isotope effect, played 
a key role in exposing their mechanism for superconduc- 
tivity. In contrast, the mechanism for superconductiv- 
ity in the cuprate is still elusive, but is believed to be 
of magnetic origin 0- Verifying this belief would re- 
quire an experiment similar to the isotope effect, namely, 
a measurement of Tc versus the magnetic interaction 
strength J, with no other structural changes in the com- 
pounds under investigation. Here we present such an ex- 
periment using the (Caa;Lai_.2:)(Bai.75_a;Lao.25+2:)Cu30j; 
(CLBLCO) system with its four different superconduct- 
ing families, for which maximum Tc (T™°^) varies by 
28%. This is a large change compared to Sn, which has 
the strongest isotope effect in nature where Tc varies only 
by 4%. For each family, we measured the Neel Temper- 
ature Tjv and the anisotropies of the magnetic interac- 
tions. This allows us to obtain the Heisenberg coupling 
J. In addition, we determine the spin glass tempera- 
ture Tg of underdoped samples. J, Tg and Tc allow us 
to generate a unified phase diagram for magnetism and 
superconductivity from no doping to over doping. We 
combine this result with a previous determination of the 
superconducting carrier density rig 0) and demonstrate 
experimentally a magnetic analog of the isotope effect. 

CLBLCO is a high temperature superconductor 
(HTSC) system that belongs to the YBasCuaOy (YBCO) 
group. Each value of x in the range of 0.1 < x < 0.4 is 
a family in the system. All families are tetragonal. The 
difference in the unit cell parameters a and c/3 between 
the two extreme families (x = 0.1 and 0.4) is 1% 0. 
Therefore, variations in T"^"'^ due to variation in ionic ra- 
dios are not relavent since they are expected to be on the 
order of 1 K, and with opposite sign to observation 0. 
Moreover, the level of disorder as detected by Ca NMR 
and Cu NQR is identical for the different families. 
The range of accessible j/'s in CLBLCO is so large that 
samples from the heavily underdoped antiferromagnetic 
(AFM) parent compounds, to the non-superconducting 
extreme overdoped, are obtainable. When varying x, 



nearly concentric domes of Tc versus oxygen doping y 
are formed, with the maximum value of Tc varying from 
58 K at X = 0.1 to 80 K at a; = 0.4 0, as demonstrated 
by the open symbols in Fig. ^ The magnetic properties 
of these compounds are determined using the zero field 
muon spin relaxation (/iSR) technique. 
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FIG. 1: Phase diagram for the 

(Caa;Lai_3;)(Bai.75_iLao.25+x)Cu30y system. 
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Fig. 121 shows typical muon polarization P{t) curves, 
at different temperatures, for three samples from the 
X = 0.1 family. At high temperatures the polarization 
curves from all samples are typical of magnetic fields em- 
anating from nuclear magnetic moments. In this case 
the time dependence of the polarization exhibits, as ex- 
pected, a Gaussian decay. As the temperature is lowered 
the sample enters a magnetic frozen phase and the polar- 
ization relaxes much more rapidly. While the transition 
from paramagnetic to the frozen state looks identical for 
all samples, the behavior at very low T is different and 
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indicates the nature of the ground state. Fig. Ola) is 
an example of an antiferromagnetic ground state. When 
the temperature decreases, long range magnetic order is 
established at ^ 377 K reflected by spontaneous oscilla- 
tions of the muon polarization. Fig. Etc) is an example 
for a spin glass (SG) transition at ~ 17 K. In this case 
the ground state consists of magnetic islands with ran- 
domly frozen electronic moments 0], and consequently, 
the polarization shows only rapid relaxation. When the 
transition is to a Neel or spin glass state, the critical tem- 
peratures are named T^r and T^, respectively. Figl^^b) 
presents an intermediate case where the sample appears 
to have two transitions. The first one starts below 240 K, 
where the fast decay in the polarization appears. Be- 
tween 160 K and 40 K there is hardly any change in the 
polarization decay, and at 30 K there is another transi- 
tion manifested in a faster decaying polarization. This 
behavior was observed in all the samples on the border 
between antiferromagnet and spin-glass in the phase di- 
agram. 
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FIG. 2: Time evolution of the muon polarization for 
(Ca2,Lai_2:)(Bai.75-a:Lao.25+a:)Cu30y from the x = 0.1 fam- 
ily close to the magnetic critical temperature, (a) A sample 
with an antiferromagnetic transition, (b) A sample with both 
an antiferromagnetic and a spin-glass transition, (c) A sam- 
ple with a spin-glass transition. The solid lines are a fit to 
Eq.lEll. 
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FIG. 3: The magnetic volume fraction extracted from the 
muon depolarization as a function of temperature for the 3 
samples shown in Fig. |5| 



In order to determine the magnetic critical transition 
temperatures, the data were fitted to a sum of two func- 
tions: a Gaussian, and a rapidly relaxing function that 
describes the magnetic field due to frozen electronic mo- 
ments, namely. 



P{t) = A„ exp(- 



2j-2 



-) + 



(1) 



-I- Am |aexp(— v^Ait) + (1 — a) exp(— ■\/A2i) cos(a;i)| . 

In this function An and Am represent the amplitudes 
(i.e. the volume fraction) of the nuclear and magnetic 
parts, respectively, and Ai, A2 are the relaxation rates of 
the magnetic part. In the fit A is determined from high 
temperatures, and the sum A^ + A^ = 1 is constant at 
all temperatures. The solid lines in Fig. |21are the fits of 
Eq.(Q to the data. Fig. O shows A^ as a function of 
temperature, for the three samples in Fig. |21 Above the 
transition, where only nuclear moments contribute. Am 
is close to zero. As the temperature decreases, the frozen 
magnetic part increases and so does Am^ at the expense 
of An- For the pure AFM and SG phases, the transition 
temperature was determined as the temperature at which 
A„i is half of the saturation value. For the samples with 
two transitions, two temperatures were determined using 
the same principle. The coexistence of Neel and spin- 
glass states indicates that the transition from one kind 
of ground state to the other is a first order quantum phase 
transition. 
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The full phase diagram of the CLBLCO compound, ob- 
tained from the measurements described above, including 
the T„ measurements in the superconducting state from 
Ref. [9j, is presented in Fig. ^ From this diagram it 
is clear that, at high doping levels, CLBLCO is a pure 
superconductor. As the oxygen doping decreases con- 
siderably from optimum, these compounds can be in a 
mixed state of magnetic islands, which freeze at Tg, and 
superconductivity in between. Upon further decrease in 
doping, the system has a pure spin glass ground state. 
Even further underdoping towards the parent compounds 
leads to abrupt replacement of the glassy state with a 
long range antiferromagnetic order with T^r ~ 400 K, 
which seems to saturate. The inset of Fig. ^ shows that 
y ^ 6.5 is a crossing point, and at higher doping, families 
with high X values have lower T/v and Tg. In contrast, at 
doping levels below y = 6.5, families with higher x values 
(and higher T^™'^^) have the higher Tn- This is the first 
indication of correlation between T™°^ and magnetic in- 
teractions. 
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FIG. 4: Upper inset: the phase diagram after the scaling de- 
scribed in the text. Lower inset: The effective anisotropies as 
a function of AP for the different families, see text for details. 
The solid lines are guides to the eye. Main figure: phase di- 
agram after both scaling and extraction of the contributions 
from anisotropies as described in the text. T^"*" = J for the 
parent compounds. 



In order to untangle the phase diagram, we use the 
scaling relation introduced in Ref. 0] that caused both 
Tc and Tg data of many cuprate families to collapse 
into a single curve 0|. This means plotting, for each 
family, Tc divided by T™"^ of the family, as a func- 
tion of AP where AP — K{x)Ay, Ay = y — yo is the 
chemical doping measured from optimum yo, and K{x) 



is the scaling parameter determined for each family so 
that all the superconducting domes collapse onto a sin- 
gle curve 0. We use K = 0.77,0.67,0.54,0.47, and 
2/0 = 7.135,7.15,7.155,7.15 for the x=0.1 to 0.4 respec- 
tivly. The results are presented in the inset of Fig. 0] We 
find that, with the exception of the x = 0.1 family in the 
AFM region, the scaling relation holds perfectly well for 
the entire phase diagram. 

Nevertheless, we would like to check whether the scal- 
ing relation can be extended to the a; = 0.1 family. We 
suspected that in x = 0.1 anisotropies might be different 
from the other families. It is well established that a pure 
2D AFM orders magnetically only at T = 0, and that T^ 
for 3D AFM is finite. Intermediate cases are described by 
anisotropic interactions where J and a J are the in and 
out of plane coupling respectively. For small a the reduc- 
tion of the magnetic order parameter M with increasing 
T is fast so that at a = the 2D limit is recovered. A 
plot the normalized order parameter a = M/Mq, where 
Mq is the order parameter at T ^ 0, versus T/T^ should 
connect (1,0) to (0,1) (see Fig.EJ in a way that is deter- 
mined only by a. At the same time J oc Tn ln(a) (up 
to In of In of a corrections) so that knowledge of a 
would lead to J. 

To test our hypothesis we measure P(t) with high tim- 
ing resolution at T — > 0, and extract a. A typical time 
dependent polarization is depicted in the inset of Fig. 
The best fit of the polarization is achieved with the func- 
tion 
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P(<) = ^Ai exp(-A,t) cos{LUit) (2) 

i=l 

with LU3 — 0, which is shown in the inset by the solid line. 
The reason for multiple frequencies is that the muons 
stop at different sites in the unit cell. Since the muon 
oscillation frequency is proportional to the local magnetic 
field it experiences, uj(T)/u!o is equivalent to a. The order 
parameter extracted from the high frequency, around few 
tens of MHz {luq 27 MHz in our case) is known to agree 
with neutron scattering determination of cr jl2l |. The 
lower frequency is believed to emerge from metastable 
muon sites and is not used for further analysis. 

In Fig. 13 we present a for the two different CLBLCO 
samples with x = 0.1 and 0.3 having the same AP that 
is marked in the inset of Fig.0J The data sets starting at 
tr = 1 is from the high frequency and will be used for com- 
parison with theory. Clearly the reduction of the mag- 
netization with increasing temperatures is not the same 
for these two samples, and therefore their anisotropies 
are different. Since a is less sensitive to increasing T in 
the x — O.l family than in the x = 0.3 family the a of 
X — 0.1 must be larger. Consequently T^ of the x — 0.1 
"turns out" to be too high due to a and not J. This 
could explain the deviation from scaling of the a; = 0.1 
family. 
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To account for the anisotropics quantitatively we as- 
sume that CLBLCO could be considered as a 2D magnet 
with weak anisotropics since the chain layers are partially 
full with oxygen even for the parent compounds. There- 
fore we use the Hamiltonian 

H = J S^ ■ Sj+5|| +axy ^ S-S-^s^^ +a± ^ • Si+s^ ) 

Z,(5|| Z,(5|| ^,(5_L 

(3) 

where 5\\ and S± are the in and out of plane neigh- 
bor spacings, respectively. We apply the self-j:onsistent 
Schwinger-boson mean- field (SBMF) theory to cal- 
culate a{aeff, t) where t = T / J , aeff — z^yU^y + z±a±, 
and the z^s are the number of neighbors. The calculation 
is done by solving simultaneously for every aeff and t 
two equations: a self consistency equation 

h = 2aeff[l~2KiA,h,t)] (4) 

and a constraint equation ensuring one Schwinger-boson 
per site 



angle of the optimally doped samples changes by 1.3 de- 
grees from X = 0.1 to a; = 0.4. Similarly the lattice con- 
stant a changes by 0.03 According to Ref. p^ . 
where J is plotted as a function of the Cu-O-Cu buckling 
angle and bond length for a variety of materials, this king 
of changes can produce a 30% variation in J. 

Determining ijv experimentally for the AFM com- 
pounds allows us to present a modified phase diagram us- 
ing r^°''. The modified phase diagram showing Tff^ ,Tg 
and Tc normalized by T^'^^ versus AP is presented in 
Fig. ^ Note, Tg is not corrected since it is determined 
by interactions between islands, and anisotropics are not 
expected to play an important role in this case. The 
large error bars for is due to the pure determination 
of In (aeff)- Nevertheless, the scaling now works for all 
four families, with T^"'^ variation of 28%, in most of the 
AFM region as well. This demonstrates that an energy 
scale Jf , unique for each family but varying between 
families, controls both T/v, Tg and Tc- 



K{A,h,t) + K{A,0,t) = 1. 



(5) 



In these equations 



^ (l + A + h) 



K{A,h,t) = 2.32 / , , 

^ ' Jo ^(A,/i,7) 

the density of states is given by 



n(w(A,/i,7),t) + l/2]p(7)d7, 
(6) 



(7) 



or its approximation m, w(A,/i,7) = 2.32((1 + A + 
/i)2_ ^2)1/2^ and n{uj,t) = [exp(w/i) - 1]"\ Finally, 
a{aeff,t) — h{aeff, t)lh(pieff, 0). tj^ and T/v are defined 
by cr(ae//, ijv) = and T^r — JIn- The theoretical order 
parameter for several aeff as a function of T /Tn is also 
presented in Fig. by the solid lines. 

When comparing theory and experiment we focus on 
the low temperature data, up to 200 K, where the the- 
ory is most accurate. We determine the best aeff and 
present them in the lower inset of Fig. 01 The solid lines 
are guides to the eye. Using the corresponding In {aeff) 
we obtain the corrected T/v defined as 



Ticor 
N 



TN/tN{aeff) 



(8) 



where T^''' — J for the parent compounds. At AP = 
—0.4 where T^r of all samples is saturated and they can 
be compared we find J = 950(60), 1130(70), 1260(80), 
and 1330(80) K for a; = 0.1, 0.2, 0.3, and 0.4 respectivly. 

This variation of J is probably due to an increasing 
buckling angle and decreasing lattice constants with in- 
creasing X. Neutron diffraction show that the buckling 
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FIG. 5: The normalized staggered magnetization as a function 
of the normalized temperature. The symbols are the experi- 
mental results, taken by measuring the oscillation frequency 
of the polarization curves. The solid lines are the theoretical 
curves plotted according to Eqs. 0]to|7| Inset: an example of 
muon polarization for a sample in the Neel state at T ^ 0, 
the solid line is a fit to Eq.lj^J 



In a previous work it was found that CLBLCO obeys 
the Uemura relation Tc oc tIj (16] , where Ug is the super- 
conducting carrier density, in both under- and overdoped 
regions 3] . Combining this fact with our present finding 
that a magnetic energy scale controls T™'*'^, suggests the 
formula 

Te = cJfns{AP) (9) 
where the constant c = 0.7(1) (using the definition 
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ns(0) = 0.08 carrier per Cu), AP = at optimal doping 
for all families, and is family independent. 

Equation @ is a magnetic equivalent of the isotope ef- 
fect Tc = cM~^/'^. We demonstrate this magnetic effect 
for 15% variation of T^'^^ without any theory. This is 
far greater than the variation due to the strongest 
isotope effect in metallic superconductors. Using the 
SBMF theory we extend the effect further to 28% of 
j-max Therefore, our data is a strong support to the 
belief that magnons are responsible for producing pair- 
ing in the cuprates. 
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